arXiv: 1502.06965vl [astro-ph.EP] 24 Feb 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 26 February 2015 (MN DTeX style file v2.2) 

Calibrating Gyrochronology using Kepler Asteroseismic targets 


Ruth Angus, 1 * Suzanne Aigrain , 1 Daniel Foreman-Mackey , 2 and Amy McQuillan 3 

1 Department of Physics, University of Oxford, UK 

"Centre for Cosmology and Particle Physics, New York University, New York, NY, USA 

3 School of Physics and Astronomy, Raymond and Beverly Sackler, Faculty of Exact Sciences, Tel Aviv University, 69978, Tel Aviv, Israel 


Draft version 2014 November 24 


ABSTRACT 

Among the available methods for dating stars, gyrochronology is a powerful one because 
it requires knowledge of only the star’s mass and rotation period. Gyrochronology relations 
have previously been calibrated using young clusters, with the Sun providing the only age 
dependence, and are therefore poorly calibrated at late ages. We used rotation period mea¬ 
surements of 310 Kepler stars with asteroseismic ages, 50 stars from the Hyades and Coma 
Berenices clusters and 6 field stars (including the Sun) with precise age measurements to cal¬ 
ibrate the gyrochronology relation, whilst fully accounting for measurement uncertainties in 
all observable quantities. We calibrated a relation of the form P = A n x {B — V — c) b , where 
P is rotation period in days, A is age in Myr, B and V are magnitudes and a, b and n are the 
free parameters of our model. We found a = 0.40+q'q 5 , b = 0.3 Aq'J^ and n = 0.55+o’g9- 
Markov Chain Monte Carlo methods were used to explore the posterior probability distribu¬ 
tion functions of the gyrochronology parameters and we carefully checked the effects of leav¬ 
ing out parts of our sample, leading us to find that no single relation beween rotation period, 
colour and age can adequately describe all the subsets of our data. The Kepler asteroseismic 
stars, cluster stars and local field stars cannot all be described by the same gyrochronology 
relation. The Kepler asteroseismic stars may be subject to observational biases, however the 
clusters show unexpected deviations from the predicted behaviour, providing concerns for the 
overall reliability of gyrochronology as a dating method. 

Key words: stars: evolution, stars: fundamental parameters, methods: statistical, stars: oscil¬ 
lations (including pulsations), stars: rotation, stars: solar-type 


1 INTRODUCTION 

1.1 Dating methods for field stars 

Many fields of astronomy rely on precise age measurements of 
Main Sequence (MS) stars. Unfortunately, age is a notoriously 
difficult quantity to measure for these stars, as observable prop¬ 
erties evolve slowly on the MS. Even with high precision spec¬ 
troscopic measurements, ages often cannot be determined accu¬ 
rately to within 20% ( |Soderblom|[2010j l. Some of the most pre¬ 
cise age measurements currently available are for stars in clusters 
where isochrones can be fitted to a coeval population with a range 
of masses, resulting in age measurements with uncertainties often 
as low as 10%. Isochronally derived field star ages, on the other 
hand, are much less precise than this, often having uncertainties of 
around 50% or more. Demand for age estimates of planet-hosting 
stars is high, but faint stars observed by Kepler are often expensive 
or impractical spectroscopic targets. Where high resolution spec¬ 
tra are unavailable, gyrochronology can be extrememly useful. Gy- 
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rochronology is a dating method that utilises the potentially pre¬ 
dictable rotation period evolution of low mass, MS stars. It requires 
only knowledge of the current rotation period—which is often eas¬ 
ily extracted from Kepler light curves—and mass (or appropriate 
proxy) of a star. However the current gyrochronology relations are 
entirely empirically calibrated and still need refining at large stel¬ 
lar ages. Kepler provides the perfect opportunity to calibrate gy¬ 
rochronology at late ages—it provides surface rotation periods of 
thousands of stars and new age estimates for hundreds of stars via 
asteroseismology^] This paper aims to use these new asteroseismic 
age measurements to improve the gyrochronology relations at late 
ages. 


1 Note that asteroseismic ages are only available for those Kepler stars 
which display high signal-to-noise Solar-like oscillations in their power 
spectra. The majority of stars that fall into this category are around the same 
temperature as. and slightly hotter than, the Sun 
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1.2 Gyrochronology 

Mass loss via a magnetised stellar wind causes magnetic braking of 
MS stars (Weber & Davis|19 67 ). A dynamo-driven magnetic field, 
generated at the tachocline (the interface between radiative and 
convective zones) locks the stellar wind to the surface of the star. 
The stellar wind corotates with the stellar surface out to the Alfven 
radius, at which point it decouples and angular momentum is lost 
from the star. The strength of the magnetic field at the surface, and 
therefore the rate of angular momentum loss, is related to rotation 
period ( jKawaler|1988| l. Due to this dependence, it is postulated that 
although stellar populations are born with a range of rotation peri¬ 
ods, the rapid rotators rapidly lose angular momentum and rotation 
periods converge onto a well-defined sequence. The timescale for 
convergence is expected to be around the age of the Hyades for 
solar mass stars: 650 Myrs (Ra dick et al.|19 87 Irwi n & Bouvier| 
|2009| >. After this time, rotation periods are thought to be indepen¬ 
dent of their initial values. Gyrochronology postulates that each star 
falls on a single three-dimensional plane described by mass, rota¬ 
tion period and age, i.e. given any two of these three properties, one 
can determine the third. The form of angular momentum evolution 
described above and calibrated in this article can only be applied 
to F, G and K MS stars. Fully convective M dwarfs have a dif¬ 
ferent dynamo-driven magnetic field: their rotation periods evolve 
over extremely long timescales and they often do not converge onto 
the mass-period-age plane, even after several Gyrs. Hot stars with 
effective temperatures greater than ~ 6250 K have shallow convec¬ 
tive zones—they are almost fully radiative—and, again, they have 
a different dynamo-driven magnetic field l [Kraft|1967| . These mas¬ 
sive stars retain their initial rotation period throughout their brief 
MS lifetimes and are therefore not suitable gyrochronology targets. 

The rate of rotation period decay for intermediate mass MS 
stars was first quantified by |Skumanich[{l972) , who observed that 
rotation period, lithium abundance and chromospheric activity de¬ 
cay was proportional to the square-root of age. Later, [Noyes et al.] 
( | 1984]> added a mass dependence to the period-activity-age relation 
after more massive stars were observed to spin down more slowly. 
The term ‘gyrochronology’ was coined by Barnes]( [2003[ who pro¬ 
posed an empirically motivated functional form for the relation be¬ 
tween period, colour and age, 

P = A n xa(B-V-c) b , (1) 


where P is rotation period (in days), A is age (in Myr), B and V 
are B and V band magnitudes respectively and a, 6, c and n are 
dimensionless free parameters. 

This gyrochronology relation was calibrated using open clus¬ 
ters, which are invaluable calibration tools since their ages are 
known relatively precisely and each cluster contains many stars of 
the same age which enables the period-mass relation to be tightly 
constrained. Unfortunately however, the majority of nearby clus¬ 
ters are young and until recently it was difficult to measure rotation 
periods for all but the youngest, most active stars (using ground- 
based observations). There is a significant dearth of precisely mea¬ 
sured ages for old stars and it is for this reason that the current 
gyrochronology relations are poorly calibrated at late ages. Barnes 
< |2007} > used 8 young open clusters aged between 30 and 650 Myrs 
to calibrate the dependence of rotation period on mass, and the Sun 
to calibrate the dependence on age. Best-fit values of n, a and b 
(c was fixed at 0.4) reported in Barnes (2007 1 are presented in ta- 
ble[T] Equation[Tjwas further calibrated by |Mamajek & Hillenbrand] 
(2008) using updated rotation period and age measurements of stars 

et al.|1995), Pleiades (Prosser et al.| 


in open clusters a Per (Pri 


|1995|[Krishnamurthi et al,| 1 998] , M34 (Meibom et al.|20lT) , and 
Hyades (Henry, private comm., |Radick et al.| 1987 1995 |Prosser| 
|et al,|[l995] |Paulson et al.||2004] . Once again, the Sun was used 
as an age anchor—a single data point specifying the shape of the 
period-age relation. Whereas |Bames| (2007| > fixed the position of 
c, the ‘colour discontinuity’ in equation]!] at 0.4, Mamajek & Hil¬ 
lenbrand (20081 allow it to be a free parameter in their model. The 
values of n, a and b, resulting from their fit are shown in table [I] 
In both of these studies a maximum likelihood fitting approach was 
used. This method relies on the assumption that uncertainties are 
Gaussian, which may not always be the case, and only takes obser¬ 
vational uncertainties on the dependent variable into account. As 
described in S|3] we adopt a fitting method that properly accounts 
for uncertainties on all three observed variables: colour, period and 
age. 

The data used in this article are described in S|2] our calibra¬ 
tion and model fitting process is outlined in S|3]and the results are 
presented and discussed in ]4] 


2 OBSERVATIONS 

The ages of 505 Kepler dwarfs and subgiants were published by 
|Chaplin et al.|(20l4| ). They made use of two global asteroseismic 
parameters—the average large frequency separation and the fre¬ 
quency of maximum oscillations power—to estimate stellar prop¬ 
erties, including the ages, with a grid-based approach that utilised 
several different search codes coupled to more than ten grids of 
stellar evolutionary models. 

The ages quoted in |Chaplin et al.| (2014) come from one of 
the grid-code combinations, with uncertainties reflecting the scat¬ 
ter between the different sets of results. [Chaplin et al.[(2014] used 
two different sets of effective temperatures: one was derived using 
an Infra-Red Flux Method (IRFM) calibration (Casagrande et al.] 
|20T0] |Silva Aguirre et al.|[2012j) and the other from a recalibra¬ 
tion of the SDSS griz filter KIC photometry by Pinsonneaul t et al.| 
1 2012] using Yale Rota ting Stellar Evolution Code (YREC) models 
( Demarque et al.|2004| ). We use the IRFM temperatures since they 
arc less dependent on metallicity, which is not well constrained for 
the asteroseismic sample, and their uncertainites are more conser¬ 
vative, however our analysis is relatively insensitive to this choice. 
87 stars in the asteroseismic catalogue have spectroscopic measure¬ 
ments of Teff , and [Fe/H], These precisely measured spectroscopic 
properties allowed more tightly constrained ages to be calculated 
for these 87 stars, which were incorporated where available. In or¬ 
der to produce a relation that predicts the age of a star using only 
observable properties, we chose to convert T e u to B — V for the 
asteroseismic sample using the relation of |Sekiguchi & Fukugita| 
(2000] , This conversion added an extra element of systematic un¬ 
certainty to our data since the metallicities provided for most of 
the asteroseismic stars are simply an average value for the field: 
—0.2 ± 0.3 dex (see e.g. |Silva Aguirre et al.|20lT] >. However, since 
the age uncertainties dominate this analysis, we do not expect this 
to have a significant impact on our results. 

The ages quoted in Chaplin et al. (2014) have typical uncer¬ 
tainties of 35%. These large uncertainties are the result of the fact 
that only approximate inferences can be made on the ages using the 
global asteroseismic parameters, neither of which has an explicit 
dependence on age. It will be possible, however, to derive more 
precise ages for a subset of these stars. By measuring the frequency 
of each oscillation mode individually, not just the global asteroseis¬ 
mic parameters, one can provide much tighter constraints on ages. 


© 0000 RAS, MNRAS 000 , 000-000 








































































Calibrating Gyrochronology using Kepler Astern seismic targets 3 


Table 1. Values of a, b, c and n in |Bames 1 2007 and |Mamajek & Hille‘nbran~d|l|20()8) and Angus et al. (2014). 


Parameter 


Barnes 


(20071 


Mamajek & Hillenbrand 


2008 


a 

b 

c 


n 


0.7725 ±0.011 
0.601 ± 0.024 
0.4 

0.5189 ± 0.0070 


0.407 ± 0.021 
0.325 ± 0.024 
0.495 ± 0.010 
0.566 ± 0.008 


Angus et al. (2014) 


0.40 
0.31 
0.45 
0.55 


+ 0.3 

0.05 

+ 0.3 

0.02 


+0.02 

- 0.09 


Ages derived from individual oscillation mode measurements can 
have uncertainties as small as 10% ( [Brown & Gilliland) 1994|[SlIval 
Aguirre et al. |2013[ >. However, measuring frequencies for individ¬ 
ual oscillation modes is a manual process and can only be applied in 
the highest signal-to-noise cases. [Chaplin et aT| (20 1 4[ i predict that 
around 150 of the 505 stars will be suitable for this individual oscil¬ 
lation mode treatment. We obtained precise ages for 42 stars from 
|Metcalfe et al.| ( |2014| l, modelled with the Asteroseismic Modeling 
Portal (AMP), with effective temperatures and metallicities from 
|Bruntt et al.|(2012) . Of the 42 stars in |Metcalfe et al.| ( f2014| >, we 
only incorporate the ‘simple stars’ (cool dwarfs) into our sample, 
ignoring the hotter F stars and more evolved subgiant stars as these 
are not expected to follow the simple gyrochronology relation. 

The Kepler light curves of the 505 asteroseismic targets dis¬ 
play quasi-periodic variations on timescales corresponding to the 
rotational periods of the stars—flux variations are produced by ac¬ 
tive regions on the stellar surface that rotate in and out of view. 
Rotation periods for 310 of these stars are published in |Garcia et al.[ 
( |2014j> who used a combination of an autocorrelation function and 
wavelet transform method to measure surface rotation. From the 
505 targets in the original sample, 310 rotation periods were reli¬ 
ably measured, 14 of which have precise asteroseismic ages from 
AMP modelling. All stars in the asteroseismic sample with rota¬ 
tion periods published by [McQuillan et al. (2014}, a lso appear in 
the rotation period catalogue of Garcia et al. ( |2014^ . There is ex¬ 
cellent agreement between rotation period measurements where the 
two catalogues overlap. Of the 114 stars which appear in both cata¬ 
logues, only the rotation periods of 4 were not consistent at the la- 
level and of these only 1, KIC 4931390, was inconsistent at greater 
than 2cr. The similarities between the two catalogues is further de¬ 
scribed in|Garcia et al.|(|2014]>. 


The asteroseismic sample covers a large range of ages, how¬ 
ever it does not provide good mass coverage across the entire range 
(see figures [T] and [2|. Few stars have temperatures below 6000 K 
(B — V ~ 0.55) and of the low mass stars, most of them are old 
(note that massive stars evolve rapidly and so we do not expect 
many in the sample). The exclusion of young, low-mass stars from 
the asteroseismic sample is due to the fact that these stars are more 
active and their power-spectra do not display high signal-to-noise 
acoustic oscillations. The omission of these as well as other types 
of stars that are not ideal asteroseismic targets from our sample 
should not bias our results. The mere lack of data in some regions 
of parameter space will not skew the best fitting model, however it 
is important to note that the resulting gyrochronology relation will 
not necessarily be descriptive of those stars not represented in this 
sample. We filled in some of the missing parameter space by adding 
50 stars to our sample from young clusters Coma Berenices (0.5 
Gyr), and the Hyades (0.625 Gyr) (see table [2]). Clusters younger 
than 0.5 Gyr often have large populations of rapid rotators that have 
not yet converged onto the gyrochronology plane, so no clusters 
younger than Coma Ber were included. Uncertainties on B — V 
colours associated with each cluster star were not provided in the 


catalogues from which rotation periods and ages were obtained. 
Since the uncertainty associated with each measurement plays such 
a key role in our analysis (see 0 - we assigned an uncertainty of 
±0.01 mag to each colour measurement, based on a realistic esti¬ 
mate of the typical uncertainties expected. The 1.1 and 0.588 Gyr 
open clusters, NGC 6811 and Praesepe, were originally included in 
our analysis. However we discovered that their period-colour rela¬ 
tions were different to those of the Hyades and Coma Ber, as well 
as to each other’s, and we therefore did not include them in our fi¬ 
nal analysis. A further 6 field stars with precise age measurements 
were added to the sample: 16 Cyg A and B, Alpha Cen A and B, 18 
Sco and, of course, the Sun (see table[3j. The entire set of 365 stars 
is shown in figures^and[2] Asteroseismic targets are shown in grey, 
with cluster and field stars in blue and the Sun in red. 


3 CALIBRATING THE GYROCHRONOLOGY 
RELATION 

3.1 The model 

The 310 asteroseismic stars in our sample have B—V colours con¬ 
verted from effective temperatures, photometric rotation periods, 
asteroseismic ages, and asteroseismic surface gravities. Each mea¬ 
surement of these properties is assumed to be independent with an 
associated Gaussian uncertainty. Not all of the cluster stars added 
to our sample have log g values; however, since we only use log g to 
separate the populations of subgiants and dwarfs (and we assume 
that the cluster stars are dwarfs) this should not affect our results. 
Following the treatment of |B antes ] |2007l l and|Mamajek & Hillen-| 
|brand| l |2008^ , rotation period was treated as the dependent variable 
throughout the modelling process. 

Hot stars and subgiants follow a different gyrochronology re¬ 
lation to MS dwarfs. Stars with effective temperatures above the 
Kraft-break, T e fr ~ 6250 K, l jKraft|19F7} do not have a thick con¬ 
vective envelope and cannot support a strong magnetic dynamo, 
so do not spin down appreciably during their MS lifetimes. Sub¬ 
giants spin down rapidly as they expand due to angular momentum 
conservation and thus diverge from the gyrochronological mass- 
period-age plane. The point in their evolution at which they depart, 
the ‘gyrochronological MS turn off’, is difficult to define. Classi¬ 
cally, MS turnoff is defined as the hottest point on a star’s path 
on the HR diagram (before it ascends the giant branch) but theory 
predicts that evolving stars begin the process of spinning down rel¬ 
atively slowly after leaving the classically defined MS ( |van Saders] 
& Pinsonneault 20131. For this reason we chose a very simple defi¬ 
nition of MS turnoff: we defined a log g boundary of 4.2 to mark the 
transition between dwarfs and giants. We tried a range of boundary 
values and found that 4.2 minimised subgiant contamination whilst 
maximising the cool dwarf sample. It was also necessary to use 
a mixture model to account for misclassified subgiants—without 
it, subgiant contamination significantly biased the resulting fit. We 
did not exclude hot stars and subgiants from our sample during the 
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Age (Gyr) 


Figure 1. Photometric rotation period vs age for 310 Kepler targets (grey circles) plus cluster and field stars (blue triangles). The Sun is shown as a red circle. 


Table 2. Clusters and References: (l) |Dobbie et al.|{~2009) . (2) |Collier Cameron et al.|{2009) , (3) |Perryman et al~|j I ( J98) . (4) |Radick et al.|(l987) . 


Cluster 

Age (Gyr) 

Number of stars 

Age ref 

Rotation period ref 

Coma Ber 

0.5 ±0.1 

28 

i 

2 

Hyades 

0.625 ± 0.05 

22 

3 

4 


modelling process, we modelled all three populations simultane¬ 
ously. This allowed for the fact that stars have some probability 
mass lying in all three regimes due to their large observational un¬ 
certainties. 

Hot MS stars were defined as those with B — V < 0.45, cor¬ 
responding to T e fi ss 6250 K for solar metallicity and log g. Since 
there is no dependence of rotation period on age for massive MS 
stars, their rotation periods were modelled as a normal distribu¬ 
tion with mean and standard deviation, Y and V, as free parame¬ 
ters. Subgiant rotation periods do depend on age and T e g. However, 
since the rotational properties of these stars are not interesting for 
the purposes of gyrochronology calibration, we also modelled them 
with a normal distribution with mean and standard deviation, Z and 
U, as free parameters. We used a mixture model for the remain¬ 
ing population of stars, consisting of cool dwarfs and misclassified, 
contaminating subgiants. The subgiants were treated as if their rota¬ 
tion periods had been drawn from a background normal distribution 
with mean and standard deviation, X and U, again inferred from 
the data and another parameter, Q, the probability of each star be¬ 
longing to that background population. The results of this analysis 
were not particularly sensitive to the choice of distribution for the 


hot stars and subgiants. These models were put in place so that we 
did not have to throw any data away and we could model everything 
at once. This was important because stars were classified according 
to their observed temperatures and log gs, which are noisy. Due to 
the large uncertainties on T e g and log g, each star therefore has 
some probability of being a subgiant, some of being a cool dwarf 
and some of being a hot star. By throwing away data, one could ac¬ 
cidentally throw away a misclassified star. We avoid this problem 
by modelling all stars simultaneously and taking a probabilistic ap¬ 
proach to classification. Inferences made about the parameters of 
the normal distributions used to model subgiants and hot stars were 
not of particular interest for the purposes of gyrochronology cali¬ 
bration. We decided to use simple normal distributions rather than 
more physically motivated models in order to remain as model- 
independent as possible. Figure [3] shows log g vs T e g for the aster- 
oseismic stars. Stars classified as cool dwarfs are shown in black, 
hot dwarfs in red and subgiants in blue. 

Ideally both the hot star (B — V < 0.45) and subgiant (log g< 
4.2 ) boundaries would be free parameters in our model. However, 
since these two populations were modelled with a relatively uncon¬ 
straining normal distributions, these boundary parameters would 
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+ + Asteroseismic stars 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

B-V 


Figure 2. Photometric rotation period vs B — V colour for the data described in figure [T] 


Table 3. Rotation periods and B — V colours for field stars with precise ages. Notes: [Davies et ah ra measured internal rotation periods for 16 Cyg A 
and B using asteroseismology. However, this is likely to be close to the surface rotation value. Rotation periods for and a Cen A and B were measured from 
variation in chromospheric emission lines. High-resolution spectropolarimetric observations were used to measure a rotation period for 18 Sco. The age of 16 
Cyg AB was measured with asteroseismology. 18 Sco’s age measurement was based chiefly on an asteroseismic analysis, however its rotation period was used 
as an additional constraint, so the age estimate is not entirely independent of rotation period for this star. An age for the a Cen system was estimated from 
spectroscopic observations with additional seismic constraints. 


ID age Prot B — V References 


16 Cyg A 

6.8 ± 0.4 

23.8 

0.66 ± 0.01 


Metcalfe et al. 

2012), 

Davies et al. 

(2014 

Moffett & Barnes 

(1979 


16 Cyg B 

6.8 ± 0.4 

23.2 ±“ 2 5 

0.66 ± 0.01 


Metcalfe et al. 

2012), 

Davies et al. 

(2014 

Moffett & Barnes 

(1979 


18 Sco 

3.7 ±0.2 

22.7 ± 0.5 

0.64 ± 0.01 



Li et al. 

2012 

Petit et al. ( 

20081, 

Vlermillioc 

(19861 



The Sun 

4.568 ± 0.001 

26.09 ± 0.1 

0.65 ± 0.001 

Bouvier & \ 

Vadhwa 

(20101 

Donahue et 

al. 

19 

16 

Cox & 

Pilachowski 

200C 


ot Cen A 

6 ± 1 

28.8 ± 2.5 

0.69 ± 0.01 

Bazot et a 

.12012 

, Yildn 

(20071, 

Hallam 

et al. 

19911, 

V4ermilliod ( 

986 


a Cen B 

6 ± 1 

38.7 ± 5.0 

0.90 ± 0.01 

Bazot et al. 

2012 

’ 

Yildiz 

)071, E 

umusque et al 

2012 

Mermilliod 

1986 



not be well behaved. Both would be pushed to higher and higher 
values until all stars were modelled with a normal distribution. In 
order to avoid this problem, we fixed these two boundaries. A future 
analysis could avoid the assumption that the gyrochronology rela¬ 
tion is infinitely narrow and assign it some intrinsic width, which 
would also be a free parameter. 

We postulate that there is a deterministic relationship between 
the ‘true’ rotation period of a cool MS star and its ‘true’ age and 
colour, described by equation[T](by ‘true’ we mean the value an ob¬ 
servable property would take, given extremely high signal-to-noise 
measurements). Rotation period also depends on log g since this 


property determines whether a star falls in the dwarf or subgiant 
regime. 

In what follows we use P to denote rotation period and define 
w = ( age, B—V, log g) as the vector of additional observational 
properties. Observations are denoted as P n and w n and the unob¬ 
served (latent), ‘true’ parameters as P n and w n for stars 1,..., N. 

In order to explore the posterior Probability Distribution Func¬ 
tions (PDFs) of the model parameters, 9, conditioned on a set of 
noisy observations, {P n , w„}, it was necessary to marginalise over 
the latent parameters, {P n , w„}. This is because the model param¬ 
eters, 9 are conditionally dependent on the ‘true" values of rota- 
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Figure 3. log g vs T e ff for the 310 asteroseismic stars. Hot dwarfs with T e fj > 6250 K and log g > 4.2 are red, subgiants with log g < 4.2 are blue. Only 
the black cool dwarfs with T e fj < 6250 K and log g < 4.2 are expected to follow the gyrochronology relation in equation^] 


tion period and colour, not the observed values themselves. The ob¬ 
served values are only conditionally dependent on the ‘true" values. 
In order therefore to infer the values of the model parameters using 
the observed values of rotation period and colour, it was necessary 
to marginalise over the latent parameters. Assuming all measure¬ 
ments are independent, the marginalised likelihood can be written 

1 v 

P({Pn, W n }|6>) = ]Q 
n =1 

The joint probability, on the right hand side of this equation, can be 
factorised as 

p(P n ,w„,P n ,w n \e) = ( 3 ) 

p(Pn | W „,9)p(P„ | P n )p( W„ | W„)p(w„), 

where we have utilised the fact that the observations, P n and 
w„ are conditionally independent of the model parameters, 8: they 
only depend on 8 through the latent parameters, P n and w n . The 
above integral can be written 

p({Pn, w„}|0) oc n!Li / P( w »|wu)dw„ (4) 

/ p(Pn|w„, 9)p(P n I P n )dP n , 

where we have used Bayes’ theorem: p(w n |w n ) oc 
p(w„|w„)p(w r ,,). The outer integral is the same for hot dwarfs, 
cool dwarfs and subgiants alike. In our model, the probability of 


p(P n , w n ,p„, w„\0)dP n d-w n . (2) 


the ‘true’ rotation period given the ‘true’ observed parameters and 
the model parameters, p(P„\w n , 8), is different in each regime 
because a different generative process is responsible for producing 
rotation periods. For cool dwarfs (B — V < 0.45 and log g > 
4.2): 


P(-Pn|w„,0) 


(1 — Q) 8 (P n — /s(w n )) 
+ Q (V2tt17 2 )~ 1/2 


exp 


C Pn~X) 2 \ 
2U 2 J ’ 


(5) 


where Q is the probability that a star is drawn from the population 
of misclassified subgiants and 

fe( w„) = A n x a(B - V - c) b (6) 


is the gyrochronology relation. For hot dwarfs (B — V < 0.45 and 
log g > 4.2 ) the generative process is: 

p(Pn | w„,0) = (V2 ttF 2 ) ' exp ^ Pn . 2 yl ^ ) O) 
and for subgiants, 

p(Pn | Wn,6») = (V2ttW 2 ^ 1 exp ^^ ^ ■ (8) 

We used hierarchical inference to account for observational 
uncertainties, following the method of Hogg et al.] {2010^ , also 
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used by Foreman-Mackey et al.| 2014) , |Rogers| ( |2014} , Morton & 
|Winn| ( [2014 1 and Demory < 2014 1 . We computed equation [5] up to 
an unimportant constant using a sampling approximation. The val¬ 
ues of P„ and w n with uncertainties, ap and cr w , reported in cat¬ 
alogues provide constraints on the posterior probability of those 
variables, under a choice of prior PDF, po(w n ). Ideally, these cat¬ 
alogues would provide posterior PDF samples, not just point esti¬ 
mates, which we could use directly, i.e. samples from 


p(w„|D„) = 


p(D,ijw„)po(w„ 

Po(D n ) 


(9) 


where p(D n |w n ) is the likelihood of the data, D„ (in this 
case, the set of Kepler lightcurves plus spectroscopic T e g and 
[Fe/H] measurements), given the model parameters, w n . p o(w„) 
is an uninformative prior probability distribution function po(w„ ) 
is an uninformative prior PDF, chosen by the fitter ( jChaplin et al.| 
|2014[ used a flat prior PDF in age and log g ). In the absence of 
posterior PDF sampleQwe generated our own from Gaussian dis¬ 
tributions with means, w n and standard deviations, <r w . J posterior 
samples were generated for each star (we used J = 500): 

w£ j) ~ p(w„ I w„), (10) 


and were used to evaluate p(w„|w n ) up to a normalisation con¬ 
stant. Using these samples we computed the marginalised likeli¬ 
hood for a single star as follows, 


p(Pn, Wn I 0) « ~T~ y\p{Pn | wj j) , 9) . (11) 

<J n 

3 = 1 

The argument inside this sum is given by the integral 

P(Pn\'Wn\d) = J P(Pn\'W n ,d)p(P n \ Pn)dP n . (12) 


Assuming that the period uncertainties are Gaussian with mean P n 
and variance a n 2 , this integral can be evaluated analytically for 
each population. For example, starting from equation [5] the result 
for the cool dwarfs is 


p(P n |w«>,0) 


1 -Q 

1 

? 

3 S 

to 

\J 2 7T (7 n 2 

2 cr n 2 

V J 


Q ( [Pn - X ] 2 

\/2VJU 2 AGKAj P y 2 [U 2 + a n 2 ] 


A similar result can be derived for the other populations. 

Finally, using these analytic results for the inner integral, the 
marginalised log-likelihood from equation[5]becomes 

logp({P„,w„} | 9) S3 

log 2 + Ell log [££l p{Pn I Wn ) , 9 

where Z is an irrelevant normalisation constant. We used emcee 
i jForeman-Mackey et al.|20l3) , an affine invariant, ensemble sam¬ 
pler Markov Chain Monte Carlo (MCMC) algorithm, to explore the 
posterior PDFs of the model parameters, 9. Flat prior PDFs were 
used for each parameter. Following the above method, a likelihood 
was computed as follows: 


2 Posterior PDF samples for asteroseismic parameters are now beginning 
to be published and will be made available in future publications. 


Table 4. Median values of a, b, c and n for individual clusters (see equation 

0 - 


Parameter 

Coma Berenices 

Hyades 

a 

b 

n 

0 417+ 0 - 08 
-0.07 

0 271 +0 ' 05 
' 1 _0.06 

0.542 zb 0.03 

0 312 +0 ' 04 
u,OiZ _ 0 . 06 

0.410Eq'q4 

0.599l°;“ 


• For each star, J samples were drawn from three normal dis¬ 
tributions: one in colour, one in age and one in log g, where the 
means and standard deviations of those distributions were the ob¬ 
served values and uncertainties. This step was performed just once 
and the following steps were performed for each likelihood evalu¬ 
ation. 

• For those samples that fell in the cool dwarf regime (B — V 
> 0.45 and log g > 4.2 ), model rotation periods were both cal¬ 
culated using equation[T]and assigned the value of parameter X. A 
Likelihood for the two model rotation periods were then evaluated 
using a Gaussian mixture model likelihood function. 

• For the samples that fell in the hot dwarf (B — V < 0.45 and 
log g > 4.2 ) and subgiant (log g < 4.2 ) regimes, likelihoods 
were calculated by comparing observed rotation periods with the 
model rotation periods for the two populations: Y and Z. 

• The total log-likelihood for each star was calculated as the sum 
of the log-likelihoods of each of the J samples. 

• Finally, the sum of individual star log-likelihoods provided the 
total log-likelihood. 


4 RESULTS AND DISCUSSION 

A gyrochronology relation was initially fit to the asteroseismic 
stars, the field stars and four clusters (Hyades, Coma Berenices, 
Praesepe and NGC 6811) all together. However, this resulted in ex¬ 
tremely multi-modal posteriors PDFs for a, b and n. After fitting a 
separate relation to various subsets of the data, it became evident 
that a lower value of b was preferred by NGC 6811, i.e. the slope 
of the log(period) — log(B — V — c) relation was shallower for 
NGC 6811 than for the Hyades, Coma Berenices and Praesepe. We 
therefore exluded NGC 6811 from our sample and attempted to fit 
a relation to the remaining data. Multi-modal posterior PDFs were 
|(13) still produced however, until Praesepe was also removed from the 
sample. The reason for this multi-modality is unclear, however we 
tentatively attribute it to Praesepe preffering a different value for the 
colour discontinuity, c, to the Hyades and Coma Berenices. We cal¬ 
culated the likelihood for Praesepe, plus the field stars (to provide 
the age dependence) with two different values of c: 0.45 and 0.5, 
^^finding a higher likelihood for c = 0.5. Since we do not fully under- 
) stand the cause of this variation in c, and in order to keep our model 
J simple, we chose to also exclude Praesepe from our final data set 
and fit a gyrochronology relation with c = 0.45 to the remaining 
data (asteroseismic stars, field stars, Hyades and Coma Berenices). 
We also fitted relations with c values ranging from 0.4 to 0.55 to 
this final data set, finding that the results were relatively insensitive 
to variations in this parameter (solar ages predicted from each best¬ 
fitting model were consistent within uncertainties). Individual fits 
to Hyades and Coma Berenices, plus the field stars are shown in 
figures |4(a)| and |4(bj| Median values of a, b and n for the two clus¬ 
ters with their 16th and 84th percentile uncertainties are provided 
in table [4] Note that none of these parameters are fully consistent 
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between the two clusters. The fact that each cluster seems to prefer 
a different value of a, b, n and c paints a concerning picture for this 
form of a gyrochronology relation which assumes one set of param¬ 
eters can be used to describe all stars. It is likely that the inability 
of equation [T] to fully describe the observed properties of our sam¬ 
ple of stars is due to the simplifying assumptions that go into this 
relationship. The fact that there is no dependence on metallicity, 
for example, may contribute to the deficiencies of this model. One 
might expect metallicity to have an effect on the rotational evolu¬ 
tion of a star since it impacts its internal structure. We have not at¬ 
tempted to calibrate a metallicity-dependent gyrochronology rela¬ 
tion, because there are no precise metallicity measurements for the 
majority of the Kepler stars, on which this study is chiefly based. In 
addition, we anticipate that the majority of stars for which this new 
gyrochronology relation will be most useful are Kepler stars, or 
stars targeted by missions like Kepler, again, the majority of which 
are unlikely to have precise metallicities. 


Table 5. Median values of the parameters describing the populations of 
non-gyrochronological stars. 


Parameter 

Median value 

u 

8±2 days 

V 

2-l±o 6 days 

w 

9-9tg;J days 

X 

14±4 days 

Y 

5.°±o 8 days 

Z 

16.1 'g days 

Q 

0 14+0.06 

0.05 



Age (Gyr) 

(a) Rotation period vs B — V — c for the Hyades (top) and Coma Berenices 
(bottom). 



O 



(b) Rotation period vs age for the Hyades (top) and Coma Berenices (bot¬ 
tom). 


Figure 4. Individual fits to the clusters and field stars. The Sun is the red 
point. The top figure shows rotation period vs l B — V — c’ for the Hyades 
(top) and Coma Berenices (bottom) with Solar and cluster age isochrones. 
The bottom figure shows rotation period vs age for the Hyades (top) and 
Coma Berenices (bottom) with the period- age relation for a constant B — V 
value of 0.65 (Solar B — V). 


The resulting highest probability values of a, b and n from 
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our fit to the final data set, with 16th and 84th percentile uncer¬ 
tainties are presented in table [T] The additional parameters of our 
model, describing the distributions of hot star and subgiant rota¬ 
tion periods, are presented in table [5] The posterior PDFs of these 
parameters were all unimodal. Note the value of Q, the parameter 
describing the fraction of misclassified subgiants is 0.14 !q'q®, i.e., 
based on our simple Tog g= 4.2’ definition of MS turn-off, which 
left only 21 stars classified as cool dwarfs, one or two of these are 
likely to be misclassified subgiants. Marginalised posterior PDFs 
for the three gyrochronology parameters are shown in figure [5] and 
the resulting relation between period and age for stars of Solar-like 
colour is shown in figure[6] The relation between period and colour 
for Solar-age stars is shown in figure [7] and for a range of ages in 
figures |8(a)|8(e)| Note that we plot rotation period vs 'B — V — c’, 
producing a straight line, in order to give a more intuitive under¬ 
standing of the quality of the fit to the data. In figures |6|to |8(e)| we 
plotted 100 draws from the posterior PDFs of the gyrochronology 
parameters as faint grey lines, in order to demonstrate the widths 
and bimodal natures of these distributions. 


Figure [5] shows that parameters a and n are correlated and 
their posterior PDFs are bimodal. The position of the second peak 
falls around a = 0.8, b = 0.34 and n = 0.44. The cause of this bi¬ 
modality is clear when looking at the faint grey lines representing 
draws from the parameter posterior PDFs in figure [6] The major¬ 
ity of these draws fall close to the best-fitting model, which passes 
neatly through the Sun and field stars, however a significant frac¬ 
tion fall below the line of best-fit, passing through the Kepler aster¬ 
oseismic stars which also mostly fall below the line. This result is 
reflected in |Garcia e£aL 1 2014jljvho model the AMP Kepler aster¬ 
oseismic stars from Metcalfe et al.[ ( |2014| l, without anchoring their 
relation to the Sun. They find that the model that best describes the 
Kepler asteroseismic stars underpredicts the rotation period of the 
Sun. The grey lines that fall beneath our best-fitting model in figure 
[6]are drawn from the smaller peak in the posterior PDFs of a and n 
and seem to describe the relation between rotation period and age 
for the Kepler asteroseismic stars. In other words, the bimodal pos¬ 
terior PDFs of a and n are produced by the disagreement between 
the Kepler asteroseismic stars and the Sun and field stars. One set 
of gyrochronology parameters is not capable of describing the Sun, 
plus field stars, and the Kepler asteroseismic stars simultaneously. 
There is more than one possible explanation for this result. Firstly, 
the asteroseismic ages could be systematically biased high. Sec¬ 
ondly, the rotation periods could be systematically underestimated. 
This could occur if, for example, these stars were rotating differen¬ 
tially and the dominant spotted regions on their surfaces were not 
equatorial. Thirdly, this could be a result of an observational bias 
produced by incomplete detection, sometimes known as Malmquist 
bias ( |Malmquist| 1920[ >. If there were a large spread in rotation pe¬ 
riods for a given stellar mass and age and. due to the detection 
bias brought about because shorter periods are easier to detect than 
longer periods, this broad range of rotation periods might be trun¬ 
cated at some upper cut-off. It would therefore appear as though 
Solar-colour stars were rotating too slowly for their age simply be¬ 
cause only rapidly rotating Kepler stars appeared in our sample. 
Kepler systematics hinder our ability to measure longer rotation 
periods, but in addition, more slowly rotating stars tend to be less 
active, with fewer surface features and are therefore more likely to 
be missing from our sample. Finally, it is possible that the Kepler 
asteroseismic stars follow a different spin-down relation to the Sun 
and field stars, perhaps due to having different metallicities. Un¬ 
fortunately it is not currently possible to identify the cause of the 


observed mismatch between the Kepler stars and the field stars and 
we leave this question for a future investigation. 

The fact that our fit is so heavily dominated by the Sun with 
its small uncertainties is, perhaps, a cause for concern. We know 
the rotation period and age of the Sun very precisely, however, if 
the Sun is not exactly representative of a typical star, the resulting 
best-fit gyrochronology relation will also not represent typical stars. 
The field stars in our sample are well represented by the best-fit gy¬ 
rochronology relation, which provides reassurance that the Sun is a 
typical star amongst this set. We did not attempt to tackle the prob¬ 
lem of how to appropriately treat the Sun as a single, highly-precise 
data point in a sample that also contains imprecise data. Instead, we 
leave this problem for future consideration. Despite the fact that a 
lot of weight is attributed to the Sun because of its small uncertain¬ 
ties, this new gyrochronology relation is still the most representa¬ 
tive, empirically calibrated relation between colour, rotation period 
and age for MS F, G and K stars to date. 

Our final, newly calibrated gyrochronology relation can be 
written in full as 

a re - l - 0'02 i a q a qi^O-OS 

P = A 055 - 0.09 x 0.40to;os(B -V- 0.45) 0 ' 31 — 0 02 , (15) 

with rotation period, P in days and age, A in Myr. An age can be 
calculated for a star with a rotation period and colour by inverting 
this relationship. Covariances between the gyrochronology param¬ 
eters should be taken into account whenever the above relation is 
used to calculate uncertainties on an age or rotation period predic¬ 
tion. In order to do this properly, posterior PDF samples should be 
incorporated into Monte-Carlo uncertainty calculation]^] 

In order to test the predictive power of the new gyrochronol¬ 
ogy relation, we inverted equation |T5] to compare previously mea¬ 
sured ages with new age predictions for the 6 field stars (see table 
[6]>. For comparison, gyrochronological ages for the field stars were 
also computed using the relations of |Bames ](j2007|t and|Mamajek &| 
Hillenbrand (2008). Uncertainties on ages predicted with the new 
relation were calculated using posterior PDF samples of the three 
parameters, a, b and n. 

The ages predicted by the three different relations are consis¬ 
tent within uncertainties, with the exception of Alpha Cen B. All 
three relations underpredict the age of the 16 Cyg system. The ro¬ 
tation periods of both 16 Cyg A and B used in this paper are astero¬ 
seismic measurements of the internal (not surface) rotation periods 
of the stars. If the stars’ cores are rotating much more rapidly than 
their surfaces, this could account for this age discrepancy. 

The goal of gyrochronology in general is to provide a means 
of predicting the age of a star given observations of its colour (or 
mass, or temperature), and rotation period. The discrepancies in 
period-colour relations between clusters, Kepler stars and nearby 
field stars in the above analysis does not bode well for the current, 
simple gyrochronology model. Until now it has been hoped that 
one single relation between period, mass and age could be used to 
describe all F, G and K MS stars. Even though we are unable to 
identify the cause of the discrepancy between the Kepler stars and 
the field stars, a different gyrochronology relation still seems to be 
required for each cluster. 

The ‘narrowness’ of the gyrochronology relation has hitherto 
been an unknown; do the three properties, age, mass and rotation 
period, truly lie on an infinitely narrow plane? Unfortunately we 
cannot fully answer this question here as the asteroseismic ages are 


3 Posterior samples for this project are available at 
https://github.com/RuthAngus/Gyro. 
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Figure 5. Marginalised likelihoods for the three gyrochronology parameters, a, b and n. Parameters a and n are correlated and their posterior PDFs are 
bimodal. The main peaks in the posterior PDFs of a and n correspond to a fit to the Sun and field stars. The smaller peak corresponds to a fit to the Kepler 
asteroseismic stars. This plot was made using triangle.py (Foreman-Mackey et al.|20l4) . 


Table 6. Field star ages 


taken from the literature, compared with predictions from this work (1), Mamajek & Hillenbrand f2008 > (2) and Barnes (2007) (3). 


Star 

Literature age (Gyr) 

Age 1 (Gyr) 

Age 2 (Gyr) 

Age 3 (Gyr) 

18 Sco 

3.7 ±0.2 

q y+2.5 
' —0.3 

q C+0.6 

q 7 + 0.8 

The Sun 

4.568 ±0.001 

4 - 6+3 

4 7+ 0 - 7 
^•'- 0.6 

4 8 +1 

Alpha Cen A 

6.0 ± 1 

5.0+3 

4 5+ 1 

5 ± 1 

Alpha Cen B 

6.0 ± 1 

6.0+;® 

4+ 1 

^-0.9 


16 Cyg A 

6.8 ±0.4 

4 n + 2 - 6 
^•- 0.7 

3 ± 2 

4 ± 2 

16 Cyg B 

6.8 ±0.4 

q 7+3.4 

3 ± 2 

4 ± 2 
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Age (Gyr) 


Figure 6. Rotation period vs age for cool dwarfs with colour within 0.1 of the Sun’s: 0.65, with gyrochronology relations of |Bames| {2007} , |Mamajek &| 
|HillenbrandH2008) and this work. The Sun is shown in red and the held stars, a Cen A, 18 Sco and 16 Cyg B from left to right, are shown in blue. The black 
points towards the lower left are cluster stars and those towards the upper right are Kepler asteroseismic stars. Each of the faint grey lines represents a sample 
drawn from the posterior probability distributions of a, b and n. Whilst most of these draws come from the large peak in the posterior PDF and fall through the 
Sun and held stars, some describe the period-age relation of the Kepler asteroseismic stars. These lines are drawn from the smaller peak in the posterior PDF. 


noisy and observational and intrinsic scatter are ambiguously inter¬ 
woven. However, if this sample of stars is subject to Malmquist bias 
then, by definition, there must be a broad range of rotation periods 
for stars of a given mass and age. A future study might include an 
extra parameter that describes the ‘width’ of the gyrochronological 
plane and attempt to detect an element of scatter above the noise 
level. Does age depend solely on rotation period and mass or do 
other variables influence stellar spin down, perhaps only becoming 
important after many Myrs? The gyrochronology model calibrated 
here neglects the effects of metallicity. This property is bound to 
have an effect on the angular momentum evolution of a star since it 
impacts internal stellar structure. A future study, investigating the 
impact of metallicity on rotational evolution will be essential for 
improving our ability to date stars using gyrochronology. 

The picture of gyrochronology will become clearer as the sam¬ 
ple of asteroseismic stars with individual mode analysis grows and 
their age uncertainties shrink. The best targets for asteroseismic 
studies are relatively inactive since they allow the easier detection 
of Solar-like oscillations. Inactive stars are also the best targets for 
gyrochronology as they are usually old and slowly rotating. How¬ 
ever these targets are not well suited for rotation period measure¬ 
ments which are most easily and precisely determined for active, 
rapidly rotating stars. K2, the repurposed Kepler mission, will pro¬ 


vide new targets for rotation studies; in particular, some relatively 
old clusters have been and will continue to be monitored by the 
spacecraft. The observing seasons of K2 are relatively short (~ 90 
days) and fields will only be observed once, so the maximum rota¬ 
tion periods measurable from K2 light curves will be considerably 
shorter than with Kepler. However these clusters may still be ex¬ 
tremely useful for gyrochronology. 


5 SUMMARY 

We have calibrated the relation between rotation period, B — 
V colour and age for MS stars with T e ff< 6250 K, using 
310 Kepler asteroseismic targets, supplemented with 6 field stars 
and 50 cluster stars. Unlike previous gyrochronology calibrations, 
our sample covers a large range of ages, observational uncertain¬ 
ties on all parameters were accounted for and the posterior PDFs of 
model parameters were explored using MCMC. Incorporating ob¬ 
servational uncertainties into the model fitting process was an es¬ 
sential component of our analysis since these uncertainties, partic¬ 
ularly on the asteroseismic ages, were large. Three populations: hot 
dwarfs, cool dwarfs and subgiants, were modelled simultaneously 
in order to account for potential misclassifications that might have 
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B-V- c 


Figure 7. Rotation period vs ‘B — V — c’ for dwarfs with age within lcr of the Sun’s age, 4.568 Gyr. The Sun is the red point. Each of the faint grey 
lines represents a sample drawn from the posterior probability distributions of a, b and n. Many samples fall below the solid black line marking the highest 
probability parameter values due to the bimodal posterior. 


arisen from large observational uncertainties. Posterior probability 
distributions of the gyrochronology parameters were explored us¬ 
ing MCMC and the impact of leaving out subsets of the sample 
was assessed, leading us to find that a single relation between rota¬ 
tion period, colour and age does not adequately describe the cluster 
data. For this reason, only the Hyades and Coma Berenices clusters 
were used to calibrate the model. Fitting equation[T]to the Hyades 
and Coma Berenices clusters, the Kepler asteroseismic stars, lo¬ 
cal field stars and the Sun resulted in bimodal posterior PDFs. The 
Kepler asteroseismic stars were not well described by the same gy¬ 
rochronology relation as the Sun, field stars and clusters. There are 
several potential explanations for the cause of this phenomenon. 
Some of the most likely are: the asteroseismic ages or the photo¬ 
metric rotation periods could be systematically biased; the rotation 
periods might be subject to Malmquist bias due to incomplete de¬ 
tection; and the different populations could have different physical 
properties (e.g. metallicity), not accounted for in the gyrochronol¬ 
ogy model. If Malmquist bias is responsible and there truely is a 
broad range of rotation periods for stars of a given age and colour, 
this still has negative implications for Gyrochronology which relies 
on the assumption that these three properties lie on a neat plane. 

The careful and well motivated modelling techniques used 
in this work allowed us to identify the potential shortcomings of 
the current gyrochronology model. Only by conducting MCMC 
parameter exploration and studying the resulting posterior PDFs 


were we able to see that one global gyrochronology relation can¬ 
not describe all subsets of our sample. The results are somewhat 
unsettling—no single model describes the sample as a whole, even 
when dropping some of the clusters and allowing for rather com¬ 
plex population subsets. This tells us that the current model is prob¬ 
ably not good enough, but it - or something very similar - is the 
model everyone has been using in gyrochronology studies so far. 
Previous studies did not model multiple cluster and field popula¬ 
tions simultaneously, or do the modelling in such a detailed way, 
so the problems didnt come through so clearly. In the future bet¬ 
ter physical models may be developed, better ages may be calcu¬ 
lated and more sensitive period search methods, not so susceptible 
to Malmquist bias may become available. When they do, the pre¬ 
scription we provide will enable the the period-colour-age relation 
to be modelled more sensitively. 

The code used in this project can be found at 
https://github.com/RuthAngus/Gyro. R.A. acknowledges funding 
from the Leverhulme Trust. A.M. acknowledges funding from 
the European Research Council under the EUs Seventh Frame¬ 
work Programme (FP7/(2007-2013)/ERC Grant Agreement No. 
291352). We wish to thank the anonymous referee for their valu¬ 
able and insightful feedback and Bill Chaplin for his thoughtful 
and detailed comments which enormously improved this paper. 
We would also like to thank Tsevi Maseh, Eric Mamajek, Daniel 
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Figure 8. P r ot vs B — V — c for dwarfs within 1 a of the reference age with the new gyrochronology relation and |Bames| ( |2007) , and |Mamajek & Hillenbrand| 
(2008) for comparison. Asteroseismic targets are black and cluster and field stars are red. Each of the faint grey lines represents a sample drawn from the 
posterior probability PDFs of a, b and n. Many samples fall below the best fit model due to the bimodal posterior PDFs. 
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Mortlock and David Hogg for useful insight and discussion that 
hugely contributed to this paper. 
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